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Abstract 

Both unconstrained and constrained minimax single facility loca- 
tion problems are considered in multidimensional space with Cheby- 
shev distance. A new solution approach is proposed within the frame- 
work of idempotent algebra to reduce the problems to solving linear 
vector equations and minimizing functionals defined on some idempo- 
tent semimodule. The approach offers a solution in a closed form that 
actually involves performing matrix- vector multiplications in terms of 
idempotent algebra for appropriate matrices and vectors. To illus- 
trate the solution procedures, numerical and graphical examples of 
two-dimensional problems are given. 

Key- Words: single facility location problem, Chebyshev distance, 
idempotent semifield, linear equation 

1 Introduction 

In the area of optimization, location problems pQ constitute a research do- 
main of continuing interest that goes back to the seventeenth century and 
classical works of P. Fermat, E. Torricelli, J. J. Sylvester, J. Steiner, and 
A. Weber. Many results in the domain are recognized as important contri- 
butions to various research fields including integer programming, combina- 
torial and graph optimization [U EJ [31 01 El [6] . 

Models and methods of idempotent algebra, i.e., the linear algebra over 
semirings with idempotent addition, are among the approaches developed to 
attack location problems. Besides that, new applications of idempotent al- 
gebra to real-life problems in engineering, manufacturing, information tech- 
nology, and other fields are arising and expanding [Tl [HI QUI [HI 021 03] . 
In terms of idempotent algebra, a range of problems that are nonlinear in 
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the ordinary sense, become linear, and so more tractable for analysis and 
solution. 

Idempotent algebra based approaches prove to be useful for solving cer- 
tain optimization problems, including idempotent analogues of linear pro- 
gramming problems and their extensions [13} I14j . as well as some location 
problems [HI I15|. [16] . A single facility location problem on a graph is exam- 
ined in [5J, where it is represented as a problem of minimizing a rational 
function of one variable in the idempotent algebra sense. However, the pro- 
posed solution of one-dimensional problems does not seem to be applicable 
in the multidimensional case. 

In [151 [T6] , a multidimensional constrained location problem on a graph 
is considered which has an objective function that takes the form of a max- 
imum of functions each depending only on one variable, and so is called 
max-separable. Though an efficient technique to solve the problem is devel- 
oped, this technique seems to be not suitable for problems with objective 
functions other than max-separable ones. 

An algebraic approach to multidimensional minimax single facility lo- 
cation problems with Chebyshev and rectilinear distances is proposed in 
[171 OS US ED] • The approach is based on new results in the spectral the- 
ory in idempotent algebra, including extremal properties of eigenvalues for 
irreducible matrices. The solution reduces to minimization of functionals de- 
fined on idempotent semimodules and involves evaluation of the eigenvalue 
and eigenvectors of a matrix. 

In this paper, we present a new algebraic solution to both unconstrained 
and constrained minimax location problems with Chebyshev distance. The 
proposed approach mainly exploits methods and techniques of solving linear 
vector equations developed in [211 |22"1 [23] rather than results from spectral 
theory. The approach offers a closed-form solution that involves performing 
a few usual matrix-vector operations. 

The rest of the paper is organized as follows. We begin with a brief 
overview of some concepts, definitions and notations of idempotent algebra, 
including idempotent semifields and semimodules, algebra of matrices, and 
linear vector equations. As another prerequisite, general solutions of linear 
vector equations are also outlined. Furthermore, solutions to some optimiza- 
tion problems are obtained to provide the basis for the analysis of location 
problems below. 

We consider an unconstrained location problem and present a new so- 
lution, which appear to be more direct and somewhat simpler than that in 
[171 [THl [191 [20] ■ New constrained location problems are then examined and 
explicit solutions to the problems are given. Finally, numerical examples 
of solving both unconstrained and constrained problem in two-dimensional 
space are shown to illustrate the results. 
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2 Definitions and Notations 



We start with a brief introduction to idempotent algebra so as to outline 
basic concepts, definitions and notations that together underlie results pre- 
sented in the paper. Further related details can be found in [TO], [U fTTj 

MM- 

2.1 Idempotent Semifield 

Consider a set X that is closed under addition © and multiplication © and 
has zero and identity 1. We assume (X, 0,1, ©, <S>) to be a commuta- 
tive semiring, where addition is idempotent and multiplication is invertible. 
Since the nonzero elements in X form a group under multiplication, the 
semiring is usually referred to as the idempotent semifield. 

The power notation is introduced in the ordinary way. We define X + = 
X \ {0}. For any x € X + and any integer p > 0, we have x — 1, P = 0, 
and 

x p = x v ~ x © x = x © x v ~ x , x~ p = {x~ 1 ) p . 

It is assumed that in the semifield, the integer power can be extended to 
the case of rational exponents, and so the semifield is taken to be radicable. 

From here on, as it is customary in ordinary algebra, we drop the mul- 
tiplication sign © . The power notation is meant in the sense of idempotent 
algebra. 

The idempotent addition induces a partial order < such that x < y if 
and only if x ffi y = y . From this definition it follows that 

x < x (By, y < x (By, 

and that both addition and multiplication are isotonic. 

Furthermore, we suppose that it is possible to complete the partial order 
into a linear order and consider the semifield as totally ordered. In what 
follows, the relation signs and the symbols min and max are thought of as 
referring to this linear order. 

As an example, consider the idempotent semifield of real numbers 

Rmax,+ = (R U {-oo}, -oo, 0, max, +). 

In the semifield lR m ax,+ , its null and identity elements are defined as 
= — oo and 1 = 0. For each x G R, there exists its inverse x _1 equal to 
—x in conventional arithmetic. For any x,y £ R, the power x y corresponds 
to the arithmetic product xy . The partial order induced by the idempotent 
addition coincides with the natural linear order on R. 
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2.2 Idempotent Semimodule 



Vector operations are routinely introduced based on the scalar addition and 
multiplication denned on X. Consider the Cartesian power X ra with its 
elements represented as column vectors. For any two vectors x = (ccj) and 
y = (i/i), and a scalar c £ X, vector addition and multiplication by scalars 
follow the rules 



{x © y\i — Xi © yi, 
Fig. Q] illustrates the operations in R m 



{cx}i 



CXi . 



. Addition of two vectors on 
the plane is subject to a "rectangle rule" that defines the sum as a diagonal of 
a rectangle formed by coordinate axes together with horizontal and vertical 
lines drawn through the end points of the vectors. Scalar multiplication is 
equivalent to a shift along a line directed at 45° angle to the axes. 



cx 2 




yi xi 




Xl CXl 



Figure 1: Vector addition (left) and scalar multiplication (right) in R 
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max, 



The set X n with these operations is a semimodule over the idempotent 
semifield X. 

Both vector addition and scalar multiplication are isotone operations in 
every arguments. 

A vector with all zero entries is referred to as the zero or null vector and 
denoted by . 

A vector is called regular if it has no zero entries. The set of all regular 
vectors of order n over X is Xlf . 

A vector y € X n is linearly dependent on vectors xi,... ,x m € X n , if 
there are scalars c\ , . . . , c m € X such that 

y = ciXi © • • • © c m x m . 

A vector y is collinear with x , if y = cx. 

For any vectors xi, . . . , x m £ X n , their linear span is defined as the set 



span(a;i, . . . , x m ) = I (J) c i x l 



ci, • • • ,c m e x > . 
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An example of the linear span of two vectors x\ and X2 in K^ax + i s 
given in Fig. [21 where span(a;i, a^) is a region between two thick lines going 
through the end points of the vectors at 45° angle to the axes of Cartesian 
coordinates. 





C2X2 f 


C\X\ © c 2 x 2 












C\X\ 










X\ 









Cl 





Figure 2: A linear span of two vectors in R^ ax + . 

For any nonzero column vector x = (xj) € X n , its pseudo- inverse is a 
row vector x = (x~), where = if X{ 7^ 0, and = otherwise. 
If x is any nonzero vector, then x~x = 1. 

For all regular vectors x, y G X" , the componentwise inequality x < y 
implies x~ > y . 

The distance between any two regular vectors x = (xi) and y = (j/j) is 
given by 

p(x,y) = yx®xy. 
In the semimodule IR^ax + 1 this distance becomes 

p(x,y) = max \xi - y l \, 

l<i<n 

and so coincides with the classical Chebyshev metric. 
2.3 Algebra of Matrices 

For conforming matrices A = (ojj), B = (bij) , and C = (Qj), matrix 
addition and multiplication together with multiplication by a scalar c £ X 
are performed according to the formulas 

{A © B}ij 

k 

-[cAj-jj = CCLij . 

A matrix with all zero entries is a zero matrix which is denoted by (D . 
A matrix is called regular if it has no zero rows. 
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Consider the set of square matrices X nxn . Any matrix that has all off- 
diagonal entries equal to is called diagonal. If a matrix has all entries 
above (below) the diagonal equal to , then it is triangular. 

The diagonal matrix with all diagonal entries equal to 1 is the identity- 
matrix denoted by /. 

For any matrix i / D and an integer p > 0, the power notation is 
routinely defined as 

A° = I, A p = A p ~ 1 A = AA p - 1 . 

With respect to matrix addition and multiplication, X nxn is an idem- 
potent semiring with identity. 

A matrix is reducible if it can be put in a block-triangular form by 
simultaneous permutations of rows and columns. Otherwise the matrix is 
irreducible. 

For every matrix A = (a^), its trace is given by 

n 

tr A = cijj. 
i=i 

2.4 Linear Operators and Linear Equations 

Every matrix A G x mxra defines a mapping from the semimodule X n to the 
semimodule X m . Since for any vectors x,i/£X n and scalar c € X, it holds 
that 

A(x (By) = Ax® Ay, A(cx) = cAx, 

the mapping can be considered as a linear operator. 

For given matrices A,Ce x mxn and vectors b, d G X m , a general linear 
equation in the unknown vector x G X n is written in the form 

Ax © b = Cx © d. 

Since there is no additive inverse, one cannot rearrange the equation in 
such a way that all terms involving the unknown x are brought to one side 
of the equation while those without x go to another side. 

Many practical problems reduce to solution of the following particular 
cases of the general equation 

Ax = d, Ax = x. 

By analogy with linear integral equations, the above two equations are 
respectively referred to as that of the first kind and that of the second 
kind. The last equation is also known in the literature as the homogeneous 
Bellman equation. 

Along with the equations, one can consider first- and second-kind in- 
equalities that have the form 

Ax < d, Ax < x. 
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3 Preliminary Results 



In this section, we present results from |21|, [22| [23] to be used below in the 
idempotent algebra based analysis of location problems. These results seem 
to have a less complicated form and provide better geometrical interpretation 
than similar algebraic solutions offered in the literature including [HJ [9j 

OS El 021 02 [Ml EH]. 

3.1 The Equation of the First Kind 

Given a matrix A € X mxn and a vector d € X m , the problem is to find all 
solutions x € X n of the equation 

Ax = d. (1) 

A solution to equation (pQ) is called the maximum solution if x\ > x 
for all solutions x of (P). 

We present a solution to equation ([1]) based on the analysis of the dis- 
tance between vectors in X m . 

Let a, represent column % in the matrix A for each % = 1, . . . , m . Con- 
sider a problem of evaluating the minimum distance in the sense of the 
metric p from the vector d to the linear span of these columns. 

Since we can write 

span(ai, . . . , a m ) = {Ax\x € X n }, 

the problem is to find vectors x € X n that minimize 

p(Ax, d) = (Ax)~d © d Ax. 

The next result gives the solution to the problem when both the matrix 
A and the vector d are regular. 

Lemma 1. Suppose A £ X mxn is a regular matrix, d € X m is a regular 
vector, and define 

A = ^{A{d-A)-)-d. 

Then it holds that 

min p(Ax, d) = A 

a;ex™ 

with the minimum attained at xq = A(d~ A)~~ . 

Fig. [3] presents examples of the mutual arrangement of the linear span 
of the columns in a matrix A and a vector d in K^ax + • m the case when 
A > 1 , the minimum distance to the nearest vector of the linear span is 
attained at the vector y = AA(d~ A)~ . 

As a consequence, we get the following result. 
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Figure 3: The set span(ai,a2) and the vector d in R^ ax+ when A = 1 
(left) and A > 1 (right). 

Theorem 1. A solution of equation ([1]) exists if and only if A = 1 . // 
solvable, the equation has the maximum solution given by 

x = (d-A)-. 

The general solution to equation dH) with arbitrary matrix A and vector 
d is considered in [2TJ [23] . 

3.2 Second-Kind Equations and Inequalities 

Suppose A € X nxn is a given matrix, and x € X™ is an unknown vector. 
We examine the equation 

Ax = x, (2) 

and the inequality 

Ax < x, (3) 

To solve both equation ([2]) and inequality ([3]) we propose an approach 
based on the use of a function Tr(^4) that takes each square matrix A to a 
scalar according to the definition 

n 

Tr(A) = 0trA m . 

m=l 

The function is exploited to examine whether the equation has a unique 
solution, many solutions, or no solution, and so may play the role of the 
determinant in conventional linear algebra. 

The solution involves evaluating matrices A*, A x , and A + . The first 
two are given by 

A* = I® A® ■■■ A 11 ^ 1 , A X =AA*. 
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The matrix A + is constructed as follows. Let af be column i in A x , 
and be its diagonal element. First we replace each column a* with that 
defined as 



Furthermore, the set of columns aj~ is reduced by removing those columns, 
if any, that are linearly dependent on others. Finally, the rest columns are 
put together to form the matrix A + . 

The general solutions to both equation and inequality of the second kind 
in the case of irreducible matrices are given by the following results. 

Theorem 2. Let A be an irreducible matrix, and x be the solution of 
equation ([2]) with the matrix A. 

Then the following statements hold: 

1) if Tr(^4) = 1, then x = A + v for any vector v of appropriate size; 

2) if Tr(A) t, then there is only the trivial solution x = 0. 

Fig. H] gives examples of solutions to equations ([2]) in R^ ax + for some 
particular matrices A = (01,02). In the left example, the solution set is 
depicted by a thick line drawn through the end point of the vector 02 . The 
solution on the right takes the form of a strip between thick lines going 
through a\ and ■ 



Figure 4: Solution sets for second-kind equations in IR^ax + • 

Lemma 2. Let A be an irreducible matrix, and x be the solution of in- 
equality ([3]) with the matrix A. 

Then the following statements hold: 

1) if Tr(A) < 1, then x = A*v for any vector v of appropriate size; 
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2) if Tr(A) > t, then there is only the trivial solution x = 0. 

Fig. [5] shows the solution of the inequality with the same matrix as in 
the left example in Fig. 21 The solution set forms a region that is bounded 
by the slanted thick lines. Note that in the case of the matrix in the right 
example of Fig. HI the solution sets of both equation and inequality coincide. 



^2- 


f X 




v\„. 











Figure 5: Solution set for a second-kind inequality in R^ ax + . 
Related results for the case of arbitrary matrices can be found in [221123]. 



4 Optimization Problems 

Given a matrix A G X mxn and vectors b,cG X m , consider a problem that 
is to find 

min ((Ax)~b © Ax). (4) 

£C€X™ 

Note that a particular case of the problem when c = b arises in the 
previous section in the context of the solution of first-kind equations. 
The next result offers a solution to problem 

Theorem 3. Suppose that A is a regular matrix, b and c are regular 
vectors, and assume that 



A = v / (A(c-A)-)-b. 

Then it holds that 



min ((Ax) 6©c Ax) = A, (5) 



where the minimum is attained at the vector 

x = A(c~Ay. 
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Proof. We first verify that A is a lower bound for the objective function in 
and then show that the function reaches the bound when x = A(c~A)~ . 
Take any vector x € X™ and consider 

r = (Ax)~b © Ax. 

From the last equality we have two inequalities 

r > c Ax, r > {Ax)~b. 

Right multiplication of the first inequality by x~ together with the ob- 
vious inequality xx~ > / give 

rx~ > c Axx > c~A. 

Furthermore, by pseudo-inverting both sides, we get the inequality x < 
r(c~A)~. Left multiplication by A followed by another application of 
pseudo-inversion leads to 

(Ax)- > r^iAic-A)-)-. 

Substitution into the second inequality results in 

r > r- l (A{c-A)-)-b = r^A 2 

and consequently, in the inequality r > A. 

It remains to verify that we have r = A when putting x = A(c~A)~ . 
Indeed, in this case we have 

r = (Ax)~b © c'Ax = A' 1 (A(c~ A)~y b © Ac~A(c~Ay = A © A = A, 
that concludes the proof. □ 

Now we present a useful property of the solution and consider particular 
cases when an extended solution set appears to exist. 

Corollary 1. Under the assumptions of Theorem^ the following statements 
hold: 

(i) any vector x that gives the minimum in problem (JH) satisfies the in- 
equality 

A l b <Ax< Ac; (6) 

(ii) if u is a solution of the equation Au = b, then the minimum in Q 
is attained at x = A~ l u; 

(Hi) if v is a solution of the equation Av = c, then the minimum in Q is 
attained at x = Av . 
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Proof. To verify statement (i) suppose a; is a vector that solves problem (jl]) . 
Prom the equality 

(Ax)~b<& c Ax = A 

we get two inequalities 

(Ax)'b < A, c Ax < A, 

and note that both vectors b and c are regular. 

With the same technique as above, from the first inequality we have 

{Ax)~ < {Ax)~bb" < Ab , 

and then arrive at the inequality Ax > A _1 b, which is the left part of (|6j). 
To get the right part, we take the second inequality and write 

Ax < cc Ax < Ac. 

To prove statements (ii) and (iii), we first show that A 2 > c~b. We 
have 

A(c~Ay < cc~A(c~A)~ = c 

and thus 

A 2 = (A(c-A)-)-b > c b. 

Now suppose Au = b and take x = A~ 1 u. With the above inequality, 
we get 

{Ax)~b © c Ax = Ab~b A l c b = A. 

In the same way we assume that Av = c and put x = Av . Substitution 
of x gives 

(Ax)~b © c~Ax = A b © Ac c = A. □ 
Now we give the solution to problem @ in a particular case when A = I . 
Theorem 4. Suppose that b and c are regular vectors, and assume that 

A = (c-b) 1 ' 2 . 

Then it holds that 

min (x~b © c~x) = A 

sex™ 

with the minimum attained at any vector x such that 

A < x < Ac. 

Proof. The statements of the Theorem directly follow from ([5]) and © pro- 
vided that A = I. □ 
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5 Unconstrained Location Problem 



In this section we examine a minimax single facility location problem with 
Chebyshev distance when no constraints are imposed on the feasible location 
area. 

Given m vectors r% = (ru, . . . , r n i) T € R n and constants w; £ IR, i = 
1, . . . , m, the problem is to determine the vectors x € (R n that provide 

min max x) + vus). (7) 

cegR™ l<i<m 

The problem is known as the unweighted Rawls problem with addends 
[3]. In accordance with the nomenclature of [2], it can be referred to as the 
multidimensional Chebyshev Messenger Boy Problem. 

It is not difficult to solve the problem on the plane by using geometric 
arguments [H [5] . Below we give a new algebraic solution that is based on 
representation of the problem in terms of the semifield (R max ,+ , and appli- 
cation of the results from the previous section. 

First we denote the objective function in problem ([7]) by ip(x) and write 

m 

Win) = Q)wip(ri,x). 

i=i 

Furthermore, we introduce the vectors 

p = win © • • • © w m r m , q = wir^ © • • • © w m r^. 

Writing the metric in terms of R mBX + , we have 
m 

<p{x) = (J) Wi(x~n © r^x) = x p © q x, 
i=i 

and then represent problem ((JJ) as 

min (x~p © q~x). (8) 

Application of Theorem [J] leads us to the following results which conform 
with that in [T71fT8] . 

Theorem 5. The minimum in problem ([8]) is given by 

A = {q-p) 1 ' 2 , 
with the minimum attained at any vector x such that 

A l p <x< Aq. 
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With the usual notation, we can reformulate the statement of Theorem [5] 
as follows. 

Corollary 2. Suppose that for each i = 1, . . . , n 

Pi = max(rji + w 1 ,...,r im + w m ), qi = mm(ru -w u ...,r im - w m ). 

Then the minimum in (JT]) is given by 

A = max(pi - gi, . . . ,p n - q n )/2, 

and attained at any vector x = (xj) with elements 

Pi - A < Xi < qi + A, i = l,...,n. 

An illustration of the solution in R 2 for two problems with all Wi = is 
provided in Fig. [6l where the given points are represented by thick dots. The 
solution involves drawing a minimal upright rectangle enclosing all points. 
The solution set is represented by thick line segments that go through the 
centers of the long sides in the rectangle between two slanting lines drawn 
through the vertices of the rectangle. 




Figure 6: Solutions in R 2 when all Wi = 0. 

Fig. offers an illustration to a problem with arbitrary constants Wi . 
Together with the solution of this problem (right picture) , we also give the 
solution to a corresponding problem that has all Wi set to zero (left picture) . 
To get the solution, we first replace each given point with two new points 
shown with empty circles, and then draw their related minimal rectangle. 
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P2 



P'l 




Figure 7: Solution to a problem with nonzero Wi . 



6 Constrained Location Problems 

Given a feasible set S € R n , we now consider a constrained location problem 

min max (p(Vj, x) + Wi). 
x£S l<i<m 

We suppose that the feasible location area S is defined either by equality- 
constraints as 



S = |£C 

or by inequality constraints as 



max [aij + Xj / 



Xi, i — 1, • • • , n f , 



Si 



x 



max (an + x,) < Xj, z 
i<i<« 



1, 



,n 



In the two-dimensional case, the constraints determine regions that are 
given by the intersection of half-planes with their border lines drawn at 45° 
angle to the coordinate axes on the plane. Specifically, the intersection can 
take the form of a strip that can be considered as a quite natural restriction 
for the feasible area in the Chebyshev Messenger Boy Problem. 

6.1 Algebraic Representation and Solution 

First we represent the objective function in terms of the idempotent semifield 

R maXj+ to get 

mm ip(x), (9) 
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where 

<p(x) = xp © q x. 

The feasible sets can be written as 

So = {;e|Ae = a;}, Si = {x\Ax < x}, 

Theorem 6. Suppose that A is an irreducible matrix with Tr(A) = 1, and 
assume that 

A = y/(A+(q-A+)-)-p. 

Then it holds that 

min (p(x) = A, 
xGSo 

where the minimum is attained at the vector 

x = AA + (q-A + )~. 
Proof. Since Tr(A) = 1 all solutions of the equation Ax = x take the form 

x = A + y 

for any vector y G R m , where m is the number of columns in the matrix 
A + , m < n. 

Substitution x = A + y into the objective function <p(x) turns the prob- 
lem ([9]) into an unconstrained problem of finding 

min ((A + y)~p © q~A + y). 

To solve this new problem we apply Theorem [3j First we evaluate 

A = y/(A+(q-A+)-)-p, 

and then conclude that the minimum in the problem is equal to A and 
attained at y = A(q~ A + )~ . 

Going back to the constrained problem, we finally get x = AA + (q~ A + )~ . 

□ 

Theorem 7. Suppose that A is an irreducible matrix with Tr(^4) < t, and 
assume that 

A = ^{A*(q-A*)-)-p. 

Then it holds that 

min <f>(x) = A, 
xeSi 

where the minimum is attained at the vector 

x = AA*(q-A*)~. 

Proof. It is sufficient to note that the solution of the inequality Ax < x is 
written in the form 

x = A*y 

for any vector y of appropriate size. All further arguments are the same as 
in the previous theorem. □ 
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6.2 An Example 

Consider a minimax single facility location problem with m = 2 , wi = w 2 
0, and 

-2 \ / 6 



Note that for this problem we have 

p = r 2 , q = r 1 . 

First assume that there are no constraints for the feasible location area. 
Following Theorem [5l we get 

A = (q-p) 1/2 = 4. 

Furthermore, we obtain 

x = A~ l p = Aq=^ 2 g 

Now consider the same problem under equality constraints 

max(xi,X2 — 3) = xi, 
max(xi — 5, X2 — 2) = 

Representation of these constraints in terms of K m ax,+ gives the vector 
equation 

where 

-3 



A 



-5 -2 



It is easy to see that the matrix A is irreducible and Tr(^4) = 1 = 0. To 
apply Theorem [6] we first calculate the matrices 

Since only the first column in i x has 1 = on the diagonal, we take 

Now we successively get 

q-A + = 2, A + (q-A+)~ = (ij) 
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Finally, we evaluate 



A = y/(A+(q-A+)-)-p = 10, 
and then find the solution 



x = AA + (q~A + )- = 




Suppose that the constraints take the form of the inequality 

Ax < x. 

According to Theorem [71 we calculate 



q A* = ( 2 -1 ) , A*(q-A*)- = 




Then we have 

A = y/(A*(q-A*)-)-p = 6, 
and eventually get the solution 



x = AA*(q~A*)- = 




We illustrate the above solutions in Fig. [8j where the given points are 
indicated with empty circles, whereas the location points are shown with 
thick dots. The solution for the unconstrained problem is located in the 
center of the squares, which represent isolines of the objective function. The 
thick dotes in bottom right vertices of squares correspond to the solution 
under the constraints Ax = x (left picture) and Ax < x (right picture) 
with the matrix A = (a%, 02) . 
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